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We elucidate how the color neutrality is harmed in the Polyakov Nambu-Jona Lasinio (PNJL) 
model at finite density within the adopted mean field approximation. Also we point out how usual 
assumption about the diagonal form of the Wilson loop may fail in the presence of the diquark 
condensate on several grounds. 



I. INTRODUCTION 



In recent years several aspects of the Quantum Chromodynamics (QCD) phase diagram in the fi — T plane have been 
investigated. Here /i denotes the quark chemical potential and T corresponds to the temperature. The studies are 
mainly related to the breaking and/or restoring of the global symmetries of the QCD lagrangian in the ground state. 
The high temperature-small chemical potential region of the phase diagram can be accessed by lattice calculations, 
sec for exainple[l|, 0, 0, 3 a-nd reference therein, as well as effective models, among others chiral perturbation theory 
(ChPT) quark meson couphng (QMC) 0, i and the Nambu-Jona Lasinio model(NJL) [1, [3 

. On the other 

hand the dual region of low temperature-large chemical potential is not accessible nowadays by lattice calculations 
because of the well known sign problem of 3-color QCD (for reviews see[ll|, [HI, In this case our knowledge of 

the QCD phase diagram is mainly based on phenomenological models at intermediate density and perturbative QCD 
at very high density. 

One of the most interesting features of the QCD phase diagram at high temperature and low density is the chiral 
crossover. The QCD chiral symmetry is explicitly broken by the bare quark masses and spontaneously broken by 
chiral condensates 1J|. In the vacuum, i.e. at /i = T = 0, the chiral condensate {'ipip) is expected to be of the order of 
— (250 MeV)'^. Lattice calculations show that there exists a temperature of the order of 200 MeV at which the chiral 
condensate smoothly drops down to a lower value, and at the same time chiral susceptibilities develop a peak, signaling 
the existence of a smooth crossover. The Goldstone bosons related to the chiral symmetry breaking are the pions, 
which arc massive below the phase transition temperature because of the bare quark masses. Their mass increases 
above the critical temperature so that they decouple from the low energy spectrum, and this can be interpreted as a 
further signal of the restoration of the chiral symmetry. The smooth crossover becomes a true phase transition when 
one artificially sets the bare quark masses to zero, so that the chiral symmetry becomes an exact symmetry of the 
QCD lagrangian, and above the critical temperature the chiral condensate is zero. 

Lattice calculations show that in the neighborhood of the chiral crossover a Polyakov loop crossovers occurs [TgI. 
[iTj . The Polyakov loop $ is an order parameter for deconfinement in the pure gauge theory [l^, [l^, [2^. In more 
detail, the expectation value of the Polyakov loop ($) is related to the free energy Fq of a heavy (thus non dynamical) 
quark by the relation ($) cx exp{—f]Fq). The calculations of ($) as a function of the temperature in lattice-QCD 
(LQCD) give the result that ($) = for T < Tq and ($) 7^ at T > Tq. Here Tq denotes the so called deconfinement 
temperature of the pure glue theory, Tq « 270 MeV. The LQCD results imply that for T < Tq the free energy of a 
static quark is infinite and therefore it can not be produced. On the other hand for T > Tq the expectation value of 
the Polyakov loop is no longer zero, thus the free energy of a static quark is finite and it can be produced. Beside this 
the correlation of a Polyakov loop and a Polyakov antiloop at finite distance, namely ($(0)i'(x)), is related to the 
interaction energy of a heavy quark and a heavy antiquark. At T < Tq the correlation grows linearly with the distance 
x, thus enforcing the confining feature of the low temperature region, while at T > Tq it becomes a screened Coulomb 
potential. The Polyakov loop is a gauge singlet but it is not invariant under an element of the Z3 group. Therefore in 
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the dcconfined phase the symmetry is spontaneously broken while in the confined phase it is unbroken. For this 
reason the Polyakov loop is often referred to as the order parameter of the confinement-deconfinement transition. 

In presence of dynamical quarks one can no longer regard $ as an order parameter. The very reason for this is 
that dynamical quarks break explicitly the Z3 symmetry and thus one expects ($) 7^ even in the previously called 
confined phase. Nevertheless the Polyakov loop is supposed to serve as an indicator of the phase transition, since 
LQCD results show that a Polyakov loop crossover occurs almost in coincidence with the chiral crossover at high 
temperature. The coincidence of the two crossovers suggests that chiral symmetry restoration and dcconfincmcnt are 
intimately related. 

These results are confirmed by effective model calculations. Among the numerous effective models the most cele- 
brated one is certainly the NJL model (for reviews see [111, [13, HI] ) ■ In the NJL model one replaces the gluon- mediated 
quark interaction by a 4-fermion interaction, and the chiral condensate is computed by solving a non linear gap equa- 



tion mimicking the BCS effective model of ordinary superconductors [2J|. Once the parameters of the model are 
adjusted in order to reproduce some phenomenological vacuum properties of QCD (meson masses, chiral condensate 
and pion decay constant) the model can be used to predict the existence of a chiral crossover at a temperature of the 
order of 200 MeV, in agreement with LQCD. Unfortunately the original version of the NJL model docs not incorporate 
gluons as dynamical degrees of freedom. As a consequence the interplay between the chiral and the Polyakov loop 
crossover can not be discussed within the model. This obstacle has been recently surmounted by the NJL model with 
the Polyakov loop (PNJL model in the following) [H, [ll] . 

In the PNJL model one introduces an effective potential U for the Polyakov loop and its conjugate. Its parameters 
are tuned to reproduce the LQCD results with no dynamical quarks, namely the deconfinement temperature Tq and 
the thermodynamics of the pure glue liquid for temperatures close to Tg. Moreover assuming the existence of a 
background temporal gluon field, the coupling of the Polyakov loop to the dynamical quarks is achieved via the QCD 
coupling of the background field and the quarks themselves. For more formal details see next Section. The PNJL 
model successfully reproduces several LQCD results at zero and small quark chemical potential, as well as results 
with imaginary chemical potential. 

We have demonstrated in [131 , that the Polyakov Nambu- Jona Lasinio (PNJL) model [25| only at nonzero baryon 
chemical potential may have some intrinsic problem of color neutrality. As stressed there, this issue does not exist 
at zero density, in no contradiction with the zero density results in [28[. However, once we introduce the baryon 
chemical potential, any PNJL-type model is no longer free from the color neutrality problem. This may be viewed as 
an artifact of the PNJL model, or possibly of the mean-field approximation scheme usually taken in the literature. In 
this article we wish to clarify how this problem can emerge at finite density, irrespective to the existence of diquark 
fields. 

II. FORMAL SETTINGS 

First we illustrate the problem at finite density but without the diquark condensate. To be confident, we here start 
with the NJL model for two- flavor quarks, combined with a matrix model in order to derive the PNJL type model [2^ . 
Apart from the gluon part, the partition function of the PNJL model is given as 

Z[^l,T] = jn^dL{x) j dqdqe~^o'^ d-'Sd^C.E{q,q;Ai)^ (1) 

with the Euclidean lagrangian density defined as 

CE{q. g, Ai) ^ q+{jo{dr + lAiix)) + m~ f)q-G [{qqf + {qilzTqf] . (2) 

L is the Wilson-loop matrix, L £ SU(3). In the hypothesis that the system is well-described by the mean-field 
approximation at large distance [soj , we replace the above group integral over the all space with a single integral over 
L, i.e., we take L{x) to be a constant over the whole space, which is an implicit assumption also in the PNJL model 

[H. 

L{x) ^ P^e^fo'^'irA.i.) P^^L L = e^^^/^, (3) 

where A4 = X)a=i with {(j)a} the real parameters. By introducing the auxiliary field a ^ 2Gqq, and then 

integrating out the quark fields we obtain 

Z[fi,T]=fdLf I?cr(a;)e-^'=«[^^''"(^)^'^1, (4) 

J Jp.B. 
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where the integral over L is replaced by the global one. The effective action in the presence of constant background 
is thus 

S,s[A4,a{x)-fi,T]= dT J -logDct,,y[{-ipJfi,A4]+m-a{x))Sp{x^y)], (5) 

where the covariant derivative is defined by ip[fJ., A4] = 7o(9,- + iA^ — fj.) — i'f ■ V^, and Sp{x — y) is the delta function 
periodic in imaginary time r. In the above equation we have assumed that in the ground state pion condensation 
does not take place, resulting in the vanishing of the pseudoscalar condensate (qj^Tq). This simplification is justified 
by our previous results that show that in the physical limit pion condensation is forbidden in the ground state of the 
model jsil,!!!,!!. 

We look for the homogeneous ground state via a saddle-point approximation. In this case, the approximated action 
is space-time independent and equal to /3yr2[744, cr; fi, T], with thermodynamic potential Q given by 

n{L, a;fx,T)^^- 2N,T J j^tr. In (l + L+e'^^-^)/^) - 2NtT J j^tr. In (l + Le-(^+^)/^) , (6) 



where E — y^p^ + M^ with M ^ m — a. As usual, the vacuum part of the above integral is divergent, so we need 
to put the sharp cutoff A in the momentum integral. We should note that in contrast to the /i = case, VL is not 
symmetric under the the gluonic charge conjugation L ^ . ^ 

Due to the global color SU(3) invariance of the thermodynamic function of the loop, we can always find g G SU(3) 
such that Q[L; ^,T] = il[gidiag9"^; = il[Ldiag; T] where Ldiag can be parameterized by 

idiag = diag.(e'^Se'^%e-''^i-'^^^). (7) 

Then integrating out the "phase" g using the property of the Haar measure dL = /-t(</3i, ip2)dipid'p2[dg], the partition 
function becomes 



where fj,{ipi,ip2) is the measure for the eigenvalue distribution of SU(3) that is proportional to the squared Vander- 
monde determinant [A(e''^S e"^^ e-*'^!-''^^)]^ [H. The correctly normalized measure becomes 

(sin(93i - (^2) - sin(2<^i + (^2) + sin((pi + 2(^2))^ 
M(<Pi,<P2) = ^ ■ (9) 

Here the variables are defined in the interval {ipi, (^2} G [0, 2tt]. Now we introduce the independent variables {$, $} 
by $(((01, (/J2) = ^trcLdiag and its conjugate $ = we have ^((^1, (^2) = ■^{'^ - 6$$ + 4$^ + 4$^ - 3<I>^$^). Then 
we may write down the full PNJL partition function as 



5^[$,$,a;/^,r] = /3y(l][$,i,a;/.,T] + Fgiuc[*,$]) . 

In the above equation r2[<E', $, cr; fi, T] ~ ^[L, cr; fi, T] whose explicit form is given by 

n[^^>,^,a■,^l,T] = ^ - 2NfT J ln(l + 3l'e-(^-^)/^ + ^^e'^iE^^^^ + e-3(s-M)/T) 
_2iVfT / ln(l + 3$e-(^+^)/^ + 3$e-2(^J+f )/t + g-3(B+M)/T) 



(10) 



(11) 



V^iuo[*I', ^] is the Polyakov-loop potential in the PNJL model specified below, which includes the phenomenological 
mass term proportional to $*$, and the potential inspired by the group theoretical constraint, i.e., the above-derived 
Vandermonde determinant. 

^ygiue[*, ^] = -t- bi{T) ln(l - 6$$ + 4$3 + 4$3 _ 3$2|,2~)^ ^^-^2) 



^ The partition function is symmetric only under the full charge conjugation acting on both the quark and gluon sectors, i.e., L — > L"*" 
together with — > —fi. 
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62 (i) and 64 (T) are the phenomenological parameters, and may be set as in [28|. 

In the next Section we show that the non-vanishing color density appears in the model only at ^ ^ 0. The point 
is that (f> and $ arc the global parameters in the PNJL model, and so are {ipi,(p2}. We should notice that by this 
simplification we can integrate over the fermions in the partition function. This simplification, however, inevitably 
introduces an artificial global color symmetry breaking in this model, which may be regarded as an artifact as we 
discuss below. 



A. A comment on the conjugated loop in the PNJL model. 

The expectation values for loop and conjugated loop can be determined in principle by [sot 

(13) 



= ^pnjl/ dVid^2 [(Re$)e-R°^^ cos(ImS'^) + (Im$)e-«<='5^ sinpmS*^) 



(Re$)e~^°^'' cos(ImS'<p) - (Im^je"^''^^ sinpmS*^) 



Here we have used the fact ImD, is odd with respect to the charge conjugation C acting on gluon sector, i.e., <f> ^ $. 
It is notable that ($) ($) but both stay real [s^l- Noting that the contribution to the sin(ImS'(^) term in Zpnjl 
from i+ configuration will be completely opposite to that from L, we could also write down the partition function 
which is manifestly real as 



^pnjlL 



[fi,T] = / d^id(^2e""'='^^cos(Im5y). (14) 



It is often said that ($) = ($) in the mean field approximation, but this is not true. It is the direct consequence 
of discarding the imaginary part of the action. In most of the literatures, C-odd piece Im is discarded, which is 
actually responsible for the equality between the expectation values of loop and antiloop. 



III. GLOBAL COLOR SYMMETRY BREAKING IN THE PNJL MODEL 



Now we turn to the problem of global color symmetry breaking. We adopt here the widely used simplification 
to discard the imaginary piece of 5^. In this way we should have ($) = ($) since the simplified action no longer 
has the C-odd part. This means that we should use the saddle point approximation with a constraint <!> = $. This 
can be done by taking a constraint tpi = —ip2 = f- This is similar to gauge fixing. Now we take the saddle point 
approximation, namely we look for the values of $ and a that minimize the thermodynamical potential: 



SReS^ 



S(p 



0, 



SReS^ 



5a 



0, 



{ip = <^0,(T = (To} 



(15) 



where ReS^ = f3V{n[^,^,a; ^i,T] +Vgiue[^,<^]) with $ = 



1+2 cos ip . 



the above equation defines the mean field values 



of $0 and ctq . This is the usual procedure in the PNJL model framework. Apparently this breaks the balance between 
red, green quark and blue quark since the Wilson loop L is asymmetric in color space within the gauge in which 



-ip2 = We name this configuration as L^Ail„, that is 



ril) 
■^diag 










e-"^ 
1 



(16) 



In this configuration we have observed a finite color density (q+Tgg) ^ ^ for /i ^ [23]; the price to pay to avoid 
this fact is to introduce the color chemical potential fig, whose self consistently computed value assures the vanishing 
of the color charge (see also Ref. [l^). ^ When the diquark condensate comes to play a role, {q^T^q) = ng 7^ at 
/ig = regardless of the presence/absence of the Polyakov loop. One could think that it does not bring any problem 
because of a residual discrete "gauge" freedom which remains even after a continuous gauge freedom is lost by the 
gauge fixing L — > gLdiag.?""'^- This remnant is nothing but a permutation of {r,g,b) quarks. 



^ At = 0, quark density and antiquark density are balanced so that there is no color density in the system even if = 0. This is 
consistent with the results obtained in [Ssl . 
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One must, of course, keep in mind that several configurations that give the same $o niay contribute to the total 
•^pnjlI^Oj '^Oj I^^T] via the integral given in Eq. (jlOp . Let us first introduce explicitly the elements of the permutation 
group, which is subgroup of SU(3). The foUowings are the elements of (1)<^ {r,g flip), (2)<^ {g,b flip), (3)<^ 
{b,r flip). 

10 

= ( 1 I , gbr = ( 1 I (17) 
1 

Using the above discrete transformations, we can produce another configuration (different gauge) which realizes the 

(2) 

same value of $. We name the {r,g) flipped configuration i^jagi that is 

^diig - ff.s^diagff-/ - I e''^ I . (18) 

From these -^jjag ^-i^d i^^agi '^^^ construct four more configurations as 

e^'^ \ / e-**^ 

(3) _ . . r(i) .-1 _ I n 1 n r (4) _ „ . r (2) 




idiig = 99>>Lli]U9;t =01 , L^tig = 9sbL^'l,9;b =010 

e-"fi I V e^'^ 

1 \ / 1 

4iig-5;..41UV = I ^^'^ ' 4iig=5;..4iig5r/- e^^ 

e"^ / \ e-"^ 



(19) 



All these configurations i^^)^ ... ^ 7^(6)} 

give the same In the original variables {(pi,(p2) parametrizing Ldiag 
as Ldiag = diag.(e''^i , e*'^^, e~"''i~"^^), these configurations correspond to {{(p,2tt — ip)i,{2'K ~ 1^9, (/3)2, (<p, 0)3, (27r — 
0)4, (0, <p)5, (0, 27r — 'f)^- This will become clearer in the diagonal color SU(3) basis, {A3, As}, writing Ldiag = 
gi</>3A3+j08A8 \^ ^]-^jg base, these configurations make a hexahedron in the color (3,8) plane, i.e., {(f>3,(j)s) = 

{{ip, 0)1, {—(p, 0)2, (f , ^(p)3, (— f , —^(p)4, (^, —^(p)5, (— f , ^ip)6}- It is worth to note that each configuration has 
a different color charge although all these configurations are degenerate in energy. Let us denote the color densities in 
configuration (1) generally as (n3,n8)i even though 713 = in this case. Then, noting that color density transforms 
by each discrete gauge transformation as {q^Taq) — > {q'^U^j^TaUabq), we obtain color density in each configuration 

("3W,«8W) = {(^13,^8)1, (-ri3, 718)2, (f + ^ng, ^713 - f )3, i-'-f + ^ns, -^n^ ~ ^)^, (f - ^7^8, -^ris - 

^^5, (— ^ — -^ns, -^^3 — ■^)6}- Wc notice that n'^^^ + n^--^ = n§ + n| is a gauge invariant quantity as it should 

One could argue that the average of color densities over all these configuration, i X]^=i("-3i "^s)^ = (0, 0) so that we 
would not need to worry about the color density. Of course this is not the case. The reason is that the thermodynamic 
potential of the system is not the average of those, that is i Yl^=i ^(i) where fi(j) is the thermodynamic potential for the 
configuration (03,08)ij i-e., (3Vfl(^i^ = Ke S^\^^^ ,pg-).. To see this more clearly, we write the mean-field approximated 
partition function as 




f^PNjL[1'o,'To;M,7^] = -T7ln Ve-^^"<" ^=4°° min (0(1), ^2), • • • , ^6))- (21) 



thus we have 



Of course, in current case with no diquark condensate, each potential is degenerate in energy although each has a 
different color charge. Nevertheless there is no valid argument to choose the average of as a definition of ri^j^j^. 

Furthermore, since the color symmetry is realized as a global symmetry in the PNJL model, these degenerated vacua 
are extremely sensitive to the infinitesimal external field which tends to align the vacuum to a specific direction, 
that is a famous spontaneous symmetry breakdown 0]. This can be easily illustrated as follows. If we introduce an 
infinitesimal external field ei which tries to align the vacuum to 17^^^ out of six directions, then we have the vacuum 
alignment to the 17'^' direction when the thermodynamic limit is taken as — > 00. Then after we switch off the 
external field ei — > 0, the vacuum remains in the (1) direction. That is 

f^PNjL[<I'o, c7o; T] = - hm hm ^ In ( ^ e-''^"(') ] = (22) 
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▲ 



(-l/2,^/3/2) 
(6) 



(l/2,^/3/2) 
(3) 



-1,0) / 
(2) /■■ 



(1,0) 

(1) 



(4) 

(-l/2,-^/3/2) 



(5) 

(l/2,-^/3/2) 



FIG. 1: The discrete gauge configurations of the Wilson line L{(f>3,(f)s) — e"^^^'*"'""^**'^^ in the (AsjAg) plane which result in 
the same value of $ = ^tVcL. The axis scale is normalized by (p such that $ = (1 + 2cos(p)/3. Each configuration yields 
a different color charge density; Denoting the color charge density of configuration (1), as (n3,ng), then we have (— n3,ng) 
for configuration (2), (f- + ^ - ^) for (3), {^'-f + _ _ If) for (4), - - for (5), 

^ — '^"^ , — ^) for (6). The value of ri^ + ?i| is the same for all configurations, and thus is gauge invariant. 



In this way, the global color symmetry can be broken in the adopted mean field PNJL model, and as a consequence, 
nonvanishing color charge remains in the system. Whichever gauge is chosen among six configurations, the statement 
that color charge does not vanish is gauge invariant since \/n'^ + is gauge invariant. Of course, we think that this 
global color symmetry breaking in the PNJL model is an artifact of the model, possibly of the adopted mean field 
approximation. 

One may still think that this kind of problem does not appear if we take as a dynamical variable $ itself rather 
than the gauge field {03,c/'8}- In this case one can not compute each color density, so we cannot judge whether the 
system is color neutral or not. Furthermore this does not work as soon as we take the diquark condensation into 
account. In that case, we need to treat {03,(/)8} as the dynamical variable since the quark part of potential is no 
longer function of a gauge invariant quantities, $ and $. Moreover, colored diquarks work as external fields which 
prefer some specific direction of vacuum among six configurations. When the critical temperature is approached, this 
external field vanishes rapidly, A —> as T — Tc —> 0~, but the vacuum alignment remains in the system. In this way, 
we can not avoid the problem of color symmetry breaking at finite density, /i ^ 0. 

Before going on let us make a clarification regarding the breaking of color symmetry in the PNJL model. It is 
well-known that the a local gauge symmetry can never be broken as stated in the Elitzur's theorem [Z^]. Thus the 
expectation value of an operator which is not invariant under the gauge transformation should vanish after all the 
configurations of gauge fields are properly taken into account into the partition function. Nevertheless, it is sometimes 
useful to fix a gauge and define a gauge dependent order parameter in weak coupling perturbative calculations just to 
extract gauge independent physical quantities, as done in the BCS theory of superconductivity. In the BCS theory, 
one assumes a non-vanishing electron pair condensate which leads to a physical gap in a electron spectrum. 

This is what happens in the PNJL model, as well as in the NJL model applied to color superconductive phases 
(see next Section). In the PNJL model within the adopted mean field approximation, one usually introduces a gauge 
dependent matrix- valued order parameter L (working in the Polyakov gauge) , and then he tries to look at the gauge- 
independent Polyakov-loop $ = ^tr^i which is, in a pure gauge limit, related to the free energy of a single static 
quark excitation in a gluonic heatbath. A color-symmetry breaking in a PNJL model might be understood in this way. 
But this scenario immediately brings the trouble we have encountered in this section, that the gauge-independent 
squared sum of color densities is nonvanishing at fi ^ unless either $ 7^ 1 or T = 0. After all, wc conclude that 
this color-symmetry breaking in a PNJL model is an artifact of too naive saddle-point approximation for the gauge 
dependent quantity L. In order to overcome this problem, we need to somehow construct the gauge-independent 
approximation scheme. 
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IV. COLOR SUPERCONDUCTIVITY WITHIN THE PNJL MODEL 



Let us finally discuss some difficulty to tackle the problem of color superconductivity in the PNJL model, although 
this problem is studied to some extent but with an implicit assumption which makes the problem much simpler 
[27! . [stI . [ssj . We wish here to clarify what is the problem behind these studies, as we already mentioned in [27| . 
Instead of Eq. ([2]) , we start with a matrix model combined with the NJL model described by the following Euclidean 
lagrangian, which would lead the system to be in the 2SC superconducting state [H, H^] at high density. 



Csiq, q, A4,) = g+(7o((9r + iAi) + if ■ W - i/i)q 



Gd 



be 



(23) 



^abc 



We introduced Sa and e^'^\ the antisymmetric matrix in color and flavor space respectively, i.e., [sa ] 
and [e^^'']ij = Cy . Introducing the bosonic field ^q^Cjr^s'^^s^^'^q, and integrating out fermions results in the 

expression of the partition function in terms of a global matrix model 



JdLe 



-S:,Jrfrrfx ^"'-J^°'-' +llogDct 



(c)_(f) 



(24) 



where the covariant derivative in the Euclidean space is de&ied by ip[fj., A^] = -7°^^ + i-f'di + {^l - ^^4)7". Here 
^4 = Sq=i 4>a^a is the temporal background gauge field related to the Wilson hne by L = e"^^^-^ , Sp{x — y) is again 
the delta function periodic in r. The transpose operation * should be supposed to work only in color and flavor space, 
not in the Dirac space. 

We adopt the saddle point approximation for S'off and make the following homogeneous ansatz for the supercon- 
ducting condensate, Aa{x) — ASas with A being real, obtaining a matrix model inspired PNJL action 



Spnjl[A*-A] ^ / ^^g-S,,[A,Ai;M,T] 



(25) 



where S^p = pVi}{A, A4; ^,T] whose expression will be given below, and we have introduced the condensate vector 
A = (0,0, A) in the color antitriplet space. We have taken care the fact that when integrated over L, the action 
SeS is only a function of the gauge invariant convolution A* • A, while it is not before the integration, i.e., 5^, or 
equivalently depend on the relative angle of L and A in the color space as we shall sec later. The real part of the 
thermodynamic potential in the presence of a constant background gauge A4 can be calculated by 



Ren[A,A4;/^,r] 



A^ 

Gd 



24 

E 

A=l 



dp 

(2^ 



Rc£'a(p) + 2rin||l 



(26) 



Ea{p) is chosen from the eigenvalues of the following non-hermitian 48 by 48 hamiltonian density matrix such that 
Rei^A is positive.^ The index A refers to the quasiquark degrees of freedom that can be enumerated as 3(color) x 
2 (flavor) x 4 (Dirac). 



H{p) 



■p^ili-iAi) A75£^"^£(f) 
A75£^'=^£(f) a • p + (Ai - iAi) 



(27) 



where a' = f^^^. It is important to note now that the direction of color symmetry breaking is gauge fixed to the 
third direction in color space. 

With this fact in mind, let us see what happens when we wish to integrate in the SU(3) matrix L. We again make 
use of the polar/phase decomposition of L as 



L = g ^Ldia^gg, dL = ^j.{ipi, ip2)dipidip2[dg] 



(28) 



^ The readers are refereed to [StI . |41|| for more detail of the derivation. 

* This is always possible since the eigenvalues of the Hamiltonian density appears in sets of eigenvalues, {Ea, —Ea} because of the 
Nambu-Gorkov degeneracy. 
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Ignoring the imaginary part of the PNJL model action as usual, the partition function becomes 

where g* is the complex conjugate of the phase matrix g. Here we clearly see that the phase g does not disappear from 
the integrand so that we can not trivially perform the integration over the phase. This means that the saddle-point 
solutions to ZpNjL with respect to the gauge fields {(^1,(^2} clearly depend on the relative angle g* between and 
A. This is in contrast to the case where we do not have the diquark condensate in the system. This means that we 
can no longer restrict ourselves to diagonal form of the temporal gauge when the orientation of diquark condensate 
is fixed to the third axis of the color space. Otherwise, we can not limit ourselves to the assumption A = (0, 0, A) 
when is parametrized by a diagonal form. This fact is already stressed in (27j . but it is usually assumed in the 
literatures that both ^4 is diagonal and A = (0, 0, A) [37| . 

What happens in that case? Let us assume without justification that the integration over [dg] can be done 

^PNjL(A^,r) = /M(^i,VP2)d^irf¥'2e-^^i^^^[^'^"°"^'^l, (30) 

where A = (0,0, A). Then we work the saddle point approximation in Lpi and Lp2. Again restricting <f> to be real, 
letting $ = 1+2 cos ^j^g partition function becomes 

^pnjl(m,I^)= E e-^^t^<="['^'<)"^'^l+^'=-[*'*l]. (31) 

{03,08}; 

A^l^j^ is the temporal background gauge field in the configuration {03, 08 }i- We have exponentiated the Vandermonde 
determinant and introduced a phenomenological Polyakov loop potential as before. In this case, the degeneracy in six 
discrete gauge configurations is lifted away. We have two degenerated minima in this case, that is configuration (1) 
and (2). This is a kind of vacuum alignment, because only the configuration (1) and (2) do not bring any additional 
symmetry breaking in ReJl but the original one, i.e., SU(3)c — > SU(2)c caused by diquark condensate A = (0,0, A). 
In the thermodynamic limit V 00, the state should be selected to be either in configuration (1) or (2) again not to 
their average. Even though both states in this case turn out to have the same color charge density {0,ns) because of 
vanishing 71,3, the physical state cannot be the statistical (symmetric) average of these two states because the PNJL 
model has only a global color symmetry, and the quantum field theory has an infinite number of degrees of freedom 
in the limit V ^ 00. This is in contrast to quantum mechanics where we can construct the symmetric states by 
average of two states. When temperature is increased, the 2SC state eventually undergoes a 2nd order transition to 
the normal quark matter phase. A decreases to zero when the critical temperature is approached, but the alignment 
of A4 in the color space remains in (1) or (2) configuration, both of which yield nonvanishing color density ng ^ 
even in the normal quark matter phase at /i 7^ 0. We thus conclude that the PNJL model without color chemical 
potentials at finite fi =^ is not free from appearance of color charge density even if diquark fields are absent. 



V. CONCLUSIONS 



We have discussed how color neutrality should be implemented in the adopted mean field PNJL model at finite 
density. We have given detailed formal arguments that explain why we need to introduce a color chemical potential at 
/Lt 7^ in order to assure color neutrality of the system. We have compared our results with those obtained in Ref. [28t . 
showing that there is not discrepancy at = (which is the regime investigated numerically in the aforementioned 
paper) . 

Moreover we have pointed out that the mean field solution can not be computed as an average over the degenerate 
ground states: as a matter of fact, spontaneous symmetry breaking means that one can choose only one of the 
degenerate ground states. In fact, the various ground states of the PNJL model can not communicate among them 
because they construct different Hilbert/Fock spaces, i.e., the overlap between those spaces vanishes after V ^ 00; 
we should choose only one of them for making calculations. This procedure, which is common both to discrete 
and to continuous symmetry, is very popular in statistical mechanics, when one for example studies effective models 
of spontaneous magnetization of Weiss domains: Nature chooses only one direction for magnetization, even if the 
ground state possesses an infinity of degenerate ground states corresponding to the infinite possible orientation of the 
magnetization. 

We have also discussed a subtle point related to the simultaneous diagonalization of the Polyakov loop L and the 
color superconductive gap matrix A. We have shown that in principle it is not correct to take both of them in a 
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diagonal form, since one needs two color rotations in order to achieve this: the color rotation that makes L diagonal 
does not necessarily make A diagonal and vice versa. Thus in principle one should work either with a diagonal L or 
with a diagonal A. Furthermore we have noticed that, defining a non superconductive state as a limit for A ^ 
of a superconductive one, the vacuum alignment makes both mathematically and physically impossible to define the 
ground state as a mean of the several degenerate ground states of the model. As a consequence, we can not avoid 
non-vanishing color charge density even in the non superconducting phase without introducing appropriate color 
chemical potentials, at least in the mean field approximation adopted here. 
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